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We briefly review the effects of selective solvation of ions in aqueous mixtures, where the ion 
densities and the composition fluctuations are strongly coupled. We then examine the surface 
tension 7 of a liquid-liquid interface in the presence of ions. We show that y can be decreased 
drastically due to the electrostatic and solvation interactions near the interface. We calculate 
how the free energy is changed due to small surface undulations in the presence of an electric 
double layer. A surface instability occurs for negative 7, which can easily be realized for 
antagonistic ion pairs near the solvent criticality. Three-dimensional simulation shows how 
the surface instability is induced. 
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1. Introduction 

In usual electrolyte theories, ions interact via the Coulombic potential in a fluid with 
a homogeneous dielectric constant s. However, in most of such theories, the microscopic 
molecular interactions between ions and solvent molecules are not explicitly considered. 
Around a microscopic ion such as Na+ or CI" in a polar fluid, the ion-dipole interaction gives 
rise to a solvation (hydration) shell composed of a number of solvent molecules (those of 
the more polar component for a mixture). The resultant solvation free energy per ion will be 
called the solvation chemical potential and will be written as //^^^j where / represents the ion 
species. It is important that fx\^^ strongly depends on the solvent density and/or the composi- 
tion for binary mixtures, with its typical values much larger than the thermal energy ksT.ln 
the supercritical region of pure water the density-dependence of ^u^^j is very large. For binary 
mixtures,^' let us introduce its derivative with respect to the composition 0: 

gi = -(ksTr'diJ^Jd(/>, (1) 

which represents the degree of selective or preferential solvation of ions. Hereafter, (p is the 
water composition in water-oil. Then > for hydrophilic ions and < for hydrophobic 
ions. For strongly hydrophilic ions, l^l^^i(<p) should decrease very steeply in a concentration 
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range < <p < (p^^^, where hydration shells are formed.'*^ Here (p'^^^ ~ 0.002 for Br~ in water- 
nitrobenzene. In this paper, we discuss the selective solvation effects in the range (p > 
assuming the presence of well-defined hydration shells even in the water-poor region, where 
the variation of fJ-[gi((f>) is milder but g,- is still large (being of order 10 for monovalent ions)."^^ 
Furthermore, for poly electrolytes and colloids, dissociation of ions can occur at their ioniz- 
able groups and the degree of dissociation sensitively depends on the local composition.'^^ 
Therefore, the ion densities and the solvent composition are strongly coupled ubiquitously 
in soft matter. Around colloids, polymers, proteins, and so on, very small heterogeneities in 
the solvent composition drastically alter the ion distributions and the local electric poten- 
tial. As our recent findings, this selective solvation induces (i) prewetting transitions on the 
surfaces of ionizable rods (polyelectrolytes)^^ and spheres (colloids), (ii) phase separation 
(precipitation) induced by a small amount of highly selective solute even outside the solvent 
coexistence curve,^^- and (iii) mesophase formation in aqueous mixtures for an antagonistic 

13-15) 

In aqueous systems, the salt effect on the surface tension y has been examined extensively 
for an air-water interface, ^^^^^ while it has not yet been well investigated for a liquid-liquid 
(oil- water) interface. For a small amount of solute, most experimental interpretations 
have been based on the Gibbs formula^^'^^^ 7 = 7o - ksTT, where 70 is the surface tension 
without solute and F represents the adsorption of solute per unit area. However, we have 
recently shown the presence of a negative electrostatic contribution yei for a small amount of 
charged solute. The resultant generalized Gibbs formula is written as,"^'^^*^^ 

r = To - ksTY + 7ei- (2) 

Using the electric flux density D and the electric field E, we have 

Tel = -J dzD-E/^n 

= -J dzp^/2, (3) 

where the integral is in the surface normal direction (|| z). In the second Une we have used 
the Poisson equation V • Z) = 4np and the electric field expression E = -VO, where p is the 
charge density and O is the electric potential. Here, D, E, and p are supposed to vanish far 
from the interface. In this paper we also assume D = sE, where e is a dielectric constant. 
The electrostatic contribution is crucial for antagonistic salts^^^ and ionic surfactants.^^ This 
result has been derived from the calculation of the excess grand potential at a planar interface. 
In this paper, we will calculate the excess free energy for small surface deformations, which 
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is indeed expressed in terms of 7 in eq.(2). 

In this paper, we also discuss the surface instability induced by an antagonistic salt con- 
sisting of hydrophilic and hydrophobic ions.'^"^'^^"^^^ An example is sodium tetraphenylbo- 
rate NaBPh4, which dissociates into hydrophilic Na+ with gi » 1 and hydrophobic BPh^ 
with -g2 » 1. The latter anion consists of four phenyl rings bonded to an ionized boron. 
Such ion pairs behave antagonistically in aqueous mixtures in the presence of composition 
heterogeneity. First, they undergo microphase separation around a liquid-liquid interface on 
the scale of the Debye screening length. This unique ion distribution produces a large elec- 
tric double layer,^'^^^ leading to a large decrease in the surface tension^^^ in agreement with 
experiments. ^^'^'^^ From x-ray reflectivity measurements, Luo et al}^^ determined such ion 
distributions around a water-nitrobenzene(NB) interface by adding BPh^ and two species of 
hydrophilic ions. Second, they interact differently with water-rich and oil-rich composition 
fluctuations, leading to mesophases (charge density waves) near the solvent criticality. In 
accord with this prediction, Sadakane et al. added a small amount of NaBPh4 to a near- 
critical mixture of D2O and 3-methylpyridine (BMP) to find a peak at an intermediate wave 
number qm(~ 0.1 ~ k) in the intensity of small-angle neutron scattering. The peak height 
was much enhanced with formation of periodic structures, (iii) Moreover, Sadakane et al. 
observed multi-lamellar (onion) structures at small volume fractions of 3MP (in D20-rich 
solvent) far from the criticality, where BPhj and solvating 3MP form charged lamellae. 
These findings demonstrate very strong hydrophobicity of BPh4. (iv) Another interesting 
phenomenon is spontaneous emulsification (formation of small water droplets) at a water- 
NB interface. ^^'^^^ It was observed when a large pure water droplet was pushed into a cell 
containing NB and antagonistic salt (tetraalkylammonium chloride). 

The organization of this paper is as follows. In Sec. 2, we will present the background 
of selective solvation. In Sec. 3, we will explain a Ginzburg-Landau model for electrolytes 
accounting for selective solvation. The surface tension formula (2) will be derived. In Sec. 4, 
we will first calculate the free energy for a slightly deformed interface and then present sim- 
ulation results on the surface instability induced by antagonistic salt. 

2. Background of selective solvation 

Let us suppose two species of ions (/ = 1, 2) with charges Zie and Z2e (Zi > 0, Z2 < 0). 
At low ion densities, the total ion chemical potentials fii in a mixture solvent are expressed as 

fii = UbT InimA]) + Z;eO + ^^„i(0), (4) 
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where is the water composition and At is the thermal de BrogUe length (but is an irrele- 
vant constant in the isothermal condition) and O is the local electric potential. For neutral 
hydrophobic particles the electrostatic term is nonexistent. The ju, is a constant in equilib- 
rium except for the interface region. We consider a liquid-liquid interface between a polar 
(water-rich) phase a and a less polar (oil-rich) phase /3 with bulk compositions cp^ and cpfj with 
(pa > (pp. Thus /isoi(0) takes different bulk values in the two phases due to its composition 
dependence. So we define 

which is called the Gibbs transfer free energy (per ion here) in electrochemistry .'^•^^^'^ The 
bulk ion densities far from the interface are written as rija in phase a and rii^ in phase yS. 
From the charge neutrality condition in the bulk regions, we require ZiUia + T/inia = and 
ZiHip + Z2n2fi = 0. The potential O tends to constants Oq, and in the bulk two phases, 
yielding a Galvani potential difference, 

AO = 0„ - (6) 

across the interface. Here O approaches its Umits on the scale of the Debye screening lengths, 

and Kp^, away from the interface, so we assume that the system extends longer than 
in phase a and k^^ in phase /3. For air-water interfaces, there are no ions in the air region, so 
= 0, however. The continuity of //, across the interface gives 

keT InimJuip) + ZieM> - A/J^^ = 0, (7) 

where / = 1, 2. After some calculations, AO is expressed as^'^^^ 

AO = [Afil^ - Ai4f,]/e(Zy + IZ2I). (8) 

Similar potential differences also appear at Uquid-solid interfaces (electrodes). The ion 
densities in the bulk two phases (in the dilute Umit) are simply related by 

(Zi + \Z2\)kBT J- 

However, if three ion species are present, the ion partitioning between two phases is much 
more complicated.^^^ 

The Gibbs transfer free energy has been determined at present for water-nitrobenzene 
(NB)^^ and water- l,2-dichloroethane(EDC)^'^^ at room temperatures, where the dielectric 
constant of NB (~ 35) is larger than that of EDC (~ 10). In the case of water-NB, the ratio 
Ayu^^/fc^r is 13.6 for Na"^, 27.1 for Ca^"^, and 11.3 for Br~ as examples of hydrophilic ions. 
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while it is -14.4 for hydrophobic BPh~. In the case of water-EDC, it is 22.7 for Na"^ and 
17.5 for Br~, while it is -14.1 for BPh~. The amplitude \lS.fx^^p\lkBT for hydrophilic ions is 
larger for EDC than for NB and is very large for multivalent ions. For H"^ (more precisely 
hydronium ions HaO"^), A//^^ assumes positive values close to those for Na^ in these two 
mixtures. 

3. Ginzburg-Landau theory of mixture electrolytes 

3.1 Ginzburg-Landau free energy 

We consider a Ginzburg-Landau free energy F for a polar binary mixture (water-oil) 
containing a small amount of a monovalent salt (Zi = 1, Z2 = -1). The ions are dilute and 
their volume fractions are negligible. The variables (p, ni, and n2 are coarse-grained ones 
varying smoothly on the molecular scale. We also neglect the image interaction/^^ though it 
was included in our previous papers.^' See our previous analysis^^ for relative importance 
between the image interaction and the solvation interaction at a Uquid-liquid interface. The F 
is the following space integral in the cell, 

.t + klV0p + ^eA. (10) 
2 Stt 

The first term /tot depends on(p,ni, and ^2 as 



/tot = m + ksT 2 ln(n,a^) - 1 - gi<p 



(11) 



If the solvent molecular volumes of the two components take a common value vq = with 
a ~ 3 A, we may assume the Bragg-Williams form^^'^^^ for the first term, 

^=cp\ncp + (l- 0)ln(l - 0) +xcP(l - 0), (12) 

where X is the interaction parameter dependent on T. The critical value of;^' is 2 without ions. 

The Ai = h(27r/mikBTy'^ in eq.(l 1) is the thermal de Broglie wavelength of the species / with 
nil being the molecular mass. The gi and g2 are the solvation coupling constants assumed to 
be constants. The coefficient C in the gradient part of eq.(lO) is of order ksT/a and should 
be determined from the surface tension data or from the scattering data. The electric field is 
written as £ = -VO. The electric potential O satisfies the Poisson equation, 

-V • eVO = Anp. (13) 

For simplicity, the dielectric constant s is assumed to depend on (p as 

S(<p) = So + Ei(f>. (14) 
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where sq and si are positive constants. A linear composition dependence of s(<p) was observed 
by Debye and Kleboth for a mixture of nitrobenzene-2,2,4-trimethylpentane.^^^ 
The solvation terms (oc g^) follow if {^l^^i(4>) (i =1,2) depend on <p linearly as 

n[^,{(j))=Ai-kBTgi(l>, (15) 

where g, are constants independent of 0. We assume the presence of well-defined solvation 
shells in the concentarion range (p > (p^^^ for hydrophilic ions (see the sentences below eq.(l)). 
The first term A, is a constant yielding a contribution linear with respect to in /tot, so it is 
irrelevant at constant ion numbers. The difference of the solvation chemical potentials in 
two-phase coexistence in eq.(5) is given by A/i^^j(0) = ksTgiAcp, where A0 = 0„ - 0^ is the 
composition difference. From eqs.(7) and (8), the Galvani potential difference and the ion 
reduction factor are expressed in terms of gjAcp as 

AO = kBT(g,-g2)A(P/2e, (16) 

"ijsMla = «2^/«2a = exp[-(gi + g2)A(p/2]. (17) 

The equilibrium interface profiles can be calculated by requiring the homogeneity of the 
chemical potentials h = SF/dcf) and /i, = 6F/6ni. Some calculations give 

h = f- CV^if, - - ksT 2 girii, (18) 

Hi = kBTilnimA',) - gi(p] + Z,eO, (19) 

where /' = df/d<p. The ion distributions are expressed in terms of <p and O in the modified 
Poisson-Boltzmann relations,^^ 

m = n] exp[g,-0 - Z,eO/^Br]. (20) 

The coefficients are determined from the conservation of the ion numbers, <«,) = 
j drriiir) = hq, where <• • •) = j dr(- ■ ■ ) denotes the space average with V being 
the cell volume. The average hq = = {112) is a given constant density. 

3.2 Interface profiles and surface tension 

First, we give typical examples of the interface profiles for ei/eo = 4/3 and e^/siksT = 
20/n. In Fig. 1, we display the ion density ni(z) and the density difference ni(z) - «2(z) 
proportional to the charge density for hydrophilic ion pairs with g\ = 4 and g2 = 2. We use 

the free energy density in eq.(12) and set;^f = 2, 1.95, and 1.92. Note that the critical value 
of;^ is shifted as^^ = 2 + (gi + g2)^Vi)no/2 in the presence of ions in the mean-field theory. 
~^ Xc-, the electric double layer at the interface diminishes. On the other hand, the ion 
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100 



Fig. 1. Normalized ion density voni(z) (left) and normalized charge density vo(ni(z) - "2(2)) (multiplied by 
10"^) (right) with = 4 and ^2 = 2 for;^f = 2 in (a), 1.95 in (b), and 1.92 in (c). 




Fig. 2. Left: Normalized ion densities, where - 3, gi - -g2 - 10, and V()ni„ = voni/3 - 2x 10^*. For this 
hydrophilic and hydrophobic ion pair, a microphase separation forming a large electric double layer is apparent. 
Right: a^Ay/ksT and a^Ayi/fc^r as functions of vowia, where yi is the first term in eq.(24) and Ayi = yi - yo 
is very close to -Fa^. This shows that the electrostatic part y^. - y - y\ - - ^ dzeE^ lin dominates over 
Ayi - -kBTY. 

distributions for a pair of strongly hydrophilic and hydrophobic ions are very singular. In the 
left panel of Fig. 2, we show the ion densities for antagonistic ion pairs with gx = -g2 = 10, 
where X = 3 and Vo^iq, = Vo^i^ = 2 x 10""^. For this hydrophilic and hydrophobic ion pair, a 
microphase separation forming a large electric double layer is apparent. 

Next, we consider the surface tension of a liquid-liquid interface in our Ginzburg-Landau 
scheme, where all the quantities depend on z. In equilibrium we minimize the grand potential 
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Q = J dzco, where the grand potential density a> is written as 

C s 

OJ = /tot + 7r0'^ + -h<p- niHi - //2«2- (21) 
2 8;r 

Using eqs.(17) and (18) we find d{o} + pO)/Jz = 2C0'0". Hereafter 0' = and <p" = 

d^<p/dz^. Thus a>(z) tends to a common constant cc>oo as z ^ +oo and 

oj = C(p'^ - p<^ + ojoo. (22) 
The surface tension y = J Jz[£t>(z) - a>oo] is then written as^'^^^ 

j = J dz\c<p'^ - . (23) 
From eq.(20) we may also devise the following formula, 

7 = J dz[m + ^(f>''-kBTn-h(f>-C„y J ^^^h^^' ^^^^ 

where C„ = /(0„) - ksTnc - h(pa andh = [f((pa) - fi.^) ~ ^sTiUa - n^)M^(p with and 
being the bulk values ofn = ni+n2.ln accord with eq.(2), the right hand side of eq.(24) may 
be expanded for small hq as 

r = ro-W- j dz^E^ + ---, (25) 

where yo is the surface tension for no = 0. In the second term, T is the preferential adsorption 
at the interface expressed as 

r r An 1 

(26) 



J dz^n -Ha - —((f> - (pa) 



An 

where An = Ua - np. Thus the second term is the solute correction in the Gibbs formula. 
The third term is the electrostatic contribution growing for antagonistic salt^'^^^ and ionic 
surfactant. ^"^^ Thus eq.(25) is the generalized Gibbs formula including the electrostatic part. 
In the right panel of Fig.2, we examine how Ay = y - yo and Ayi = y j - yo are decreased with 

increasing riia for = -g2 = 10, where yi is the first term on the right hand side of eq.(24). 

1 /2 

We notice the following, (i) The changes Ay and Ayi are both proportional to n^^ at small 
VoHia. Here |Ay|/(vonia)^^^ is of order unity, so Ay is appreciable even for very small voUia- 
This salt-density dependence was observed for a water-air interface. ^^'^^^ (ii) The electrostatic 
part ye = - J dzsE^/Sn is known to be important in this case from comparison between Ayi 
and Ay = Ayi -I- ye. (iii) We confirm that the modified Gibbs relation Ayi = -ksTT holds 
excellently. 
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4. Interface deformations 

4.1 Free energy change 

We examine how the free energy changes for small deformations of the interface position 
z = ^{x,y). To this end, we superimpose small deviations 6(p and 6ni on the equilibrium 
interface profiles, <p = cf>(z) and n, = respectively. For simplicity, we replace (p(z) and 
by (p(z - - respectively. For small f . we set 

6<p = -(f>\zX, Sm = -n'^izX, (27) 

where cp' = dcp/dz and n'^ = dni/dz. Then the deviation of the dielectric constant 6s = -e\(f)'^ 
and that of the charge density 6p = -p'dC, are linear in ^, where p' = e{n\ - n'r^. However, the 
electric potential is the solution of the Poisson equation (12) and its deviation is expanded as 
Oi + 02 + 0{^^) up to order where Oi ~ ^ and O2 ~ To first and second orders in ^, 
eq.(13) yields 

V • [feV€) + eVOi] = Anp'^, (28) 

V • [JeVOi + eV02] = 0. (29) 

Hereafter O = <b{z) denotes the unpertubed potential for ^ = 0. The second-order change in 
the electrostatic free energy IS written as 



6Fe = f Jrf— IVOiP + — f& 
j L8;r 4n\ 



50i d^2' 
+ s- 



dz dz 

= J drs\V^i\^/Sn. (30) 

If we multiply eq.(29) by O and integrate over space, we find J dr(!>'(Ssd(!?i/dz + ed(!?2/dz) = 
0. Thus, in eq.(30), the second line follows from the first line. Furthermore, taking the deriva- 
tive of eq.(18) with respect to z gives fcp' - C4>"' - eiO'0"/4;T - ksT Y^iSiK = 0, where 
/" = d^f/dip^, (j)'" = d^(p/dz\ O' = d^/dz, and O" = d^^/dz^. We also find the rela- 
tion Y,i(6ni)^/ni = Y,igi6ni6(p + 5pQ>'^ from eq.(20). Using these relations, we calculate the 
second-order deviation of the free energy change F as 

6F = ^ cfr[^0'2|V^^P + ^Va)i • V^], (31) 

where V_l = (d/dx, d/dy) and we introduce the combination ^ = Oi -I- ^O' linear in ^. From 
eq.(28) (p is related to ^ by 

V-sV(p = sO'Vl^, (32) 
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SO ^ = for homogeneous ^. Using the surface tension expression in eq.(23), we obtain 

6F = Jjr^||V^^P + ^dr^\'^^\\ (33) 

where dr^ = dxdy in the first term is the integral in the xy plane. For homogeneous ^, we 
have 6F = since the deviations (26) represent a uniform translation of the interface. 

Let us consider the two-dimensional Fourier transform = j dr^e~'^'*''-^(r^}, where 
k = (kjc, ky) is the lateral wave vector. Then eq.(33) is rewritten as 

where A^fc^ arises from the second term in eq.(33). To calculate the Fourier transform 

(^^(z) = f drj_e~'^'''^(p{rj_,z) is introduced. From eq.(32) we may set (fij^^iz) = -k'^Gkiz)^]^, 
where G^(z) is a function of z satisfying 

(V,£V,-e^2)G, = eO'. (35) 

In terms of Gk, A^t is expressed as 

A, = -k^ J dz s(zW(z)Gk(z)/4n. (36) 

In particular, we consider the long wavelength limit with k = \k\ <^ , Ka, and /c^, where ^ is 
the interface thickness and Kq. and ac^ are the Debye wave numbers in the bulk two phases. In 
the monovalent case they are given by 

Kp = {%nnpe^lspkBTyi\ (37) 

where Sa and sp are the dielectric constants in the two phases. In this case, for \z\ larger than 
^, /c~\ and /c^\ we have Gi{z) = Woe~'^/(Sa + sp)k, where 

Wo = -J dze^'. (38) 

Thus Ajt is proportional to k as 

A, = WlklAniSa + ep). (39) 

If ^ is shorter than and k^^, the potential <I)(z) may be calculated from the nonlinear 
Poisson-Boltzmann equation,^^^ leading to 



k T 

Wo = e/jAO + — (£„ - ep) In 



1 + be^""^ 



(40) 



1 +^e-A[//2 

where b = spKp/SaKa and AU = eAO/^gr. If AU is small, we have Wq = (I + b)~^(sp + 
bsa)A^, which simply follows from the Debye-Hiickel approximation. In the simple case 
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Fig. 3. (Color online) Time-evolution of vqHi (left top), voni (right top), vo(ni - na) (left bottom), and U - 
e^lksT (right top), where a small amount of hydrophilic cations (^i = 12) and hydrophobic anions (g2 - -12) 
are added in the water-rich region z > at f = 0. An electric double layer is established on a time scale of 100?o, 
leading to a nnegative surface tension. 



Sa = Sj3, we have Wq = e^AO and Ak = kBT(AU)^igP /An, where £b = ISaksT is the 
Bjerrum length common in the two phases. 

4.2 Surface mode 

We consider the time development of the surface perturbation ^]^{t). In the stable case, 
the linear damping rate Q-u is given by Q-k = yklrj for «; rjok^/p, where rjo is the shear 
viscosity. ^"^^ Here we neglect the acceleration term in the momentum equation (as in eq.(44) 
below). In our case with ions, some calculations give 

= (r + A,)fc/77o. (41) 

This result can also be used when y < is realized by a change in the temperature or in the 
salt amoun, where the surface perturbation grows in time at long wavelengths. 

4.3 Numerical results for surface instability induced by an antagonistic salt 

As the right panel of Fig. 2 indicates, the surface tension y can be made negative with 
addition of an antagonistic salt. Negativity of y can be realized particularly near the solvent 



J. Phys. Soc. Jpn. 



DRAFT 




Fig. 4. (Color online) Slow growth of interface deformations induced by antagonistic ion pairs for f/fo = 
250, 500, 750, and 1000, where the colors represent U according to the color bar at the bottom. 



critical point even for small uq, where the solvent surface tension jq becomes very small. 
When y is slightly negative, we compare the two terms in eq.(34) to find a characteristic 
wave number in the early stage, 

K-An{s, + Sp)\y\IWl (42) 
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where we use eqs.(33) and (38). For Sa = ep, we have kc ~ 8;r^g|y|/(^Br|A?7|^). Near the 
solvent critical point, we have yo ~ 0.1A:gr/^^ and kc ~ 47r^B|7/yol/(^A0)^(gi - g2)^. We 
expect that the surface undulations grow around this wave number in the early stage. 

In this section, we numerically examine the resultant interface instability, where we ini- 
tially add a small amount of antagonistic ion pairs with gi = -g2 = 12 to the water side at 
<ni) = <R2) = no = 0.003vq^ The initial interface position is at the middle of the cell, so 
ni(z, 0) and n2(z, 0) are equal to 2no above the interface and below it at f = 0. The static 
parameters employed are;^^ = 2.1, C = ksTa'^, si = 0, and £b = e^/sokBT = 3a. The Debye 
wave number is then k = (S^too^b)^^^ = 0.48a~^ Space is measured in units of a = Vq^. 

The water composition (p and the ion densities obey^'^^ 

^+V-(0v) = -^V2/i, (43) 
ot ksT 

^ + V-(n,v) = -^V-n,V//,-. (44) 
ot ksT 

where Lq is the kinetic coefficient, Dq is the diffusion constant common to the cations and the 
anions, and h and //,■ are defined by eqs.(17) and (18). Neglecting the acceleration term, we 
determine the velocity field v using the Stokes approximation, 

T^oV^v = Vpi + + 2 mVni. (45) 

We introduce pi to ensure the incompressibility condition V-v = 0. We set rjoLo/kBT = 0.16a'* 
and Do = a^/to. Time is measured in units of 

to = a^/Lo. (46) 

The right hand side of eq.(45) is also written as V • n, where n is the stress tensor arising 
from the fluctuations of <p and n,. Here the total free energy F in eq.(lO) satisfies dF/dt < 
with these equations (if the boundary effect arising from the surface free energy is neglected). 
We integrated eq.(43) and (44) using eq.(45) at each step on a 64 x 64 x 64 lattice under the 
periodic boundary condition in the x and y directions. The system length is then L = 64a. At 
the top z = L/2 and the bottom z = -L/2, we required = -dQ>ldz ' and drii/dz = 0. The 
system was initially in a two-phase with a planar interface at z = separating the two phases. 

In Fig.3, we show the time development of ni, n2, ni - n2 in the very early stage, and 
O, which illustrate estabUshment of the electric double layer in the time region t < 100?o. In 
this stage, the profiles of (p and n, are nearly one-dimensional and the surface undulations are 
negligible. We can see establishment of a large electric double layer, where the electric field 
is of order ksT/e. The integral in eq.(23) is time-dependent here and may be regarded as an 
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effective surface tension. It is about 0.2kBT/a^ at ? = and tends to a negative value of order 
-O.OSkBT/a^ for t >50to in this case. Here, AU ~ 5 and eq.(42) gives kc ~ 0.3a~^. 

In Fig.4, we display slower time-evolution of interface deformations after the establish- 
ment of the electric double layer induced by antagonistic ion pairs at t/to = 250,500,750, 
and 1000, where the domain size is of order 2n/kc in accord with eq.(42). While the linear 
growth theory holds for t <500?(), the system is in the nonlinear growth regime with large 
surface disturbances at t = lOOO^o. Note that this final time is of the order of the diffusion 
time to = L^/4Do ~ 10^?o- For t :^ to the system should tend to a bicontinuous mesophase 
displayed in our previous paper. 

5. Summary 

We stress the crucial role of the selective solvation of a solute in phase transitions of var- 
ious soft materials, which should be relevant in understanding a wide range of mysterious 
phenomena in water. Particularly remarkable in polar binary mixtures are mesophase forma- 
tion induced by an antagonistic salt^"^^ and precipitation induced by a one-sided solute (a 
salt composed of hydrophilic cations and anions and a neutral hydrophobic solute). ^^'^^^ Re- 
garding the problem of mesophase formation, our theory is still insufficient and cannot well 
explain the complicated phase behavior disclosed by the experiments. '^"^^^ 

We propose experiments of the surface instability. For example, it should be induced 
around a droplet of pure water (oil) inserted in a bulk oil-rich (water-rich) region in the im- 
misible condition, where the droplet or the surrounding region contains a small amount of 
an antagonistic salt. The previous experiments of droplet formation^^'^^^ were performed in 
the strong segregation condition far from the solvent criticality and, as a result, the surface 
tension might have not been negative even after the formation of an electric double layer. 
Thus, experiments should also be performed on near-critical aqueous mixtures in two-phase 
coexistence, where a small amount of an antagonistic salt is added to either of the water-rich 
or oil-rich region as in Sec.4 of this paper. 
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